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SUMMARY 


The formulation given in this paper was developed to calculate 
flutter frequency and flutter speed for a problem with two degrees of 
freedom. Two different methods of evaluating the flutter determinant 
were used so that the results from each method could be compared. 
Although the method was developed for a particular problem applica- 
tion, it is sufficiently general to solve any flutter system that can 
be characterized by two degrees of freedom. 


INTRODUCTION 

The solution method (hereinafter referred to as Program A4868) 
was developed for computing flutter characteristics of a Pylon (strut) 
that was installed in the NASA Langley Research Center VSTOL tunnel. 
This paper presents the development of two methods used to calculate 
flutter frequency and flutter speed for a problem with two degrees of 
freedom. Also given in the paper are a program flow diagram, partial 
program listing, and a sample problem with input and output for the two 
different methods. A comparison of solution results between the two 
methods is also given. The flutter equations are derived and provided 
in the Appendix. 
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SYMBOLS LIST 


a 

b 


c 

c 



El 

GJ 


g 

I 


1/K 

L 

z 

M 


m 


S 


V 


'DIV 


Distance between elastic axis and center of mass, m (in.) 

Distance between elastic axis and center of pressure, 
m (in.) 

Local chord, m (in.) 

Local chord in streamwise direction, m (in.) 

Aerodynamic lift coefficient 
Aerodynamic lift coefficient slope, 1/rad. 

Aerodynamic pitching moment 

Bending stiffness, N - m 2 (lb - in 2 ) 

2 2 

Torsional stiffness, N - m (lb - in ) 

Structural damping 

Product of distributed mass moment of inertia about center 

2 2 

of mass axis and airfoil span, N - sec - m 
(lb - sec 2 - in 2 ) 

Dimensionless parameter, — 

Aerodynamic loading on differential segment, newton (lb.) 
Airfoil span, m (in.) (see sketch) 

Aerodynamic moment on differential segment, joule (in - lb) 

2 

Distributed mass, ( — ) 

in* 

2 2 

Reference area, m (in. ) 

Flow velocity, m/sec (in. /sec.) 

Divergence speed, m/sec. (in. /sec.) 

iii 



X 

z 


a 

Y 

P 
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03 
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03 . 


Flutter speed, m/sec. (in. /sec.) 

Airfoil spanwise coordinate m (in.) 

Displacement normal to airfoil section, m. (in.) 
Angle of attack, rad. 

Airfoil (strut) sweep angle, rad. 

Density of airstream, N/m^ (lb./in^) 

Phase angle, rad. 

Circular frequency, rad. /sec. 

Uncoupled torsional natural frequency, rad. /sec. 

Flutter frequency, rad. /sec. 

Uncoupled bending natural frequency, rad. /sec. 


Matrix Symbols 


[D] 2 

[M] 2 

[K] 2 

[RJ 2 

(q> 2 


x 2 aerodynamics matrix associated with displacements 
x 2 mass matrix 
x 2 stiffness matrix 
x 2 aerodynamic matrix 
x 1 matrix of the generalized coordinates 


iv 



FORMULATION AND SOLUTION OF FLUTTER DETERMINANT 


Program A4868 uses two methods to solve the flutter determinant. 
The equations for a two-degree-of- freedom system (see Appendix A) 
are given by 


[M] {q.} + (1 + i g) [K] {q^ - V 2 [D] {q.} + V [R] {q^ = {0} j = } 2 

’(1) 

Note: The variables and coefficient matrices of equation (1) are 
defined in Appendix A. 

Assuming solutions of the form 


qj = q x e 


itot 


q 2 = q 2 e 


i(wt + ((>) 


gives 

[- w 2 [M] + (1 + i g) [K] - V 2 [D] + i u V [R]] {q.} = {0} ‘ = ] g (2) 
for a nontrivial solution, 

|- a) [M] + (1 + i g) [K] - V 2 [D] + i to V [R]| = 0 + 1 0 


Noting that = D^ = D 2 j = 0, the flutter determinant becomes 


- u> M^ + K n + i g K n + 1 u V R 


11 


l ll I " “ 2 M 12 " V 2 D 12 + i w V R 12 


to M 21 + i co V R 21 


- ^22 + *^22 + 9 ^22 

- V 2 D 22 + i (0 V R 22 


= o 
(3) 
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Solution Method 1. 


The first method used to solve for the flutter condition Involves 
the solution of the flutter determinant In terms of a complex variable 
and the parameter for a given Mach No. This solutibn is developed 
as follows: 

From equation (1) we can write 


( - a> 2 Mji + (1 + i g) + i a) V q^ + 

( - io 2 M 12 - V 2 D 12 + i a) V R 12 ) q 2 = 0 


Dividing by wr M^ yields: 

K n ' * i V R;. 


Defining £ = £ and ^ = » y 2 


Equation (5) becomes: 


( - 1 + (^) (1 + i g) + i (j^) q^ + 

/ M 12 / 1\ 2 °12 . /lx R 12 x q 

’ 1 V q 2 " 0 


(4) 


(5) 


+ 


( 6 ) 
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In a similar manner the 2nd equation from (1) can be written as 
( - co 2 M^ + i to V R 21 ) qj + ( - to 2 M 22 + + i g 


V 2 D 22 + i to V R^) q 2 = 0 


(7) 


Dividing equation (7) by “ H 22 and introducing 


K 22 _ 2 1 _ V 

M 22 w a * K a) 


yields 


/ ^1 . . 1 k 21 x „ . / . “a . . x 

( "M^ +1 K q l + ( " 1 ~J~ (1 + 1 


1 R 21 


^ £| +i * S |><2 = ° 


( 8 ) 


CO 

Multiplying the 2nd term of equation (6) by (— ) 

u a 


yields 


( . 1+( ^, 2 (!♦,,)♦ 1* J|U) qi 


a 


+ ( /l^ 2 P 12 , . 1 R 12\ n 

( “ M n " 'K' M n K q 2 " 0 


(9) 


CO 

Making the substitution Z = (~) [1 + i g] in equations (8) and (9), 

and letting 



4 


. R 11 

b . H 12 
*17 ’ 

c Jjl 
M 11 

.il 

M 22 

r . ^1 

”22 ' 

H - ° 22 
M 22 


D = 


M 


12 _ 

11 


J = 


_ *22 


M 22 


0 ) 

(-*) 

“a 


-2 

0) 


the flutter determinant becomes: 


((Z w 2 - 1) + i A) 
( - E + i f ) 


( - B - (£) 2 C) + i A) 
(Z - 1 - (^) 2 H + i A) 


= 0 


( 10 ) 


Expanding the flutter determinant yields a complex polynomial Z which 
can be solved using the quadratic formula: 

2 

I 2 Z 2 - ((5 2 (1 - H (i) ) - 1) + (A (i) + Jj 2 J (|)) 1) Z 

2 2 2 2 

+ (1 + H (i) - A J (i) - E B - C E (|) + D F (£) ) + 

2 2 

(£ (- A - A H (£) - J + E D + B F + C F (£) ) i (11) 

Let X = -^ and 
A A = w 2 

B B = (to 2 ( 1 - H X 2 ) - 1) + (A X + w 2 J X) i 
CC=(1+HX 2 -AJX 2 -EB-CEX 2 +DFX 2 )+ 
(X(-A-AHX 2 -J+ED+BF+CF (X) 2 ) i 
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Then the two complex roots are 


7 d — - B B +lf (B B) - 4 A A C C BN4AACC 

2 A A ZN 2 AA 

The flutter conditions can be determined by computing values of Z for 
assumed values of 1/K. The root for which the imaginary part of Z changes 
sign gives the value of 1/K at which flutter is possible. By definition 


Re Z = (-£) , and Im Z = (- ± ) g. 


Therefore: g = ^ 


From a plot of g vs. (— ) with 1/K as a parameter, then the critical 

to 

(0 

condition (^p)p is determined for the actual value of g for the 
structure. Knowing w for the corresponding value of 1/K (at which 
curve intercepts actual g value) then the flutter speed is obtained 
from the relation 


02 ) 
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Solution Method 2. 


The second method used to find the flutter speed involves plotting 
the solutions of the real and imaginary parts of the flutter determinant. 
Rewriting equation (9) as 


(- 



1 + (— ) 
'to ' 



(i + i g) 



+ i Mj2 (h 
+ ( ‘ “ ( F 


2 D 
M 


12 

11 


+ 1 <2 = 0 


12 


(13) 


and equation (8) 


-i i Ro i w 

( ' B 22 + 1 ^ )<ll + ( ' l + ( ^ ) <1 + l9) 




(14) 


0) 


Define X = (— ) 

to 


Then the flutter determinant may be written as 


- 1 + x (i> (1 + 1 g) + i (£) flU- 


M 


12 


2 D 


11 


M 


11 


/lx" 12 . • ,1* 
" ( K } 1 ¥ 


M 


li 

11 


“21 


+ i 


1 R 21 


M 22 ' ‘ ^^22 


-1 + X(i t i g ,-(1, 2 ^ + I(i,^ 


(15) 
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Making the substitutions 


y 

w = — ■ L 
0) 

a 


E = 


A = 


% 

m 22 • 


F = 


R 11 
H u ' 

M 12 

B = -n— , 
M 11 

c = ^ 

Mn 

R 21 
h 22 • 

H - ° 22 
H -H^> 

, R 22 
- m 22 


R 

M 


12 

11 


and expanding the determinant yields: 

2 2 2 
(5 (1 - g 2 ) + 2 u g i ) X 2 + (- 1 - 5 (1 + (£) H + g J (£)) + 

2 7 2 

(- g - 5 (g + (£) H g - X J) + A (£) + A £ g) 1) X + 

2 2 2 2 
((1 + (£) H - (±) A J - E B - E C (i) + F D (1) ) + 

2 2 
(£ {- A - A H (|) - J + E D + F B + F C (£) ) ) i ) = 0 

Separating the real and imaginary parts into two equations yields: 

2 2 2 
(w (1 - g 2 )) x 2 + (-1.-5 (1 + (£) H + g J (£)) X 

2 2 2 2 

+ ((1 + (£) H - (i) A J - E B - E C - E C (£) + F D (£) ) = 0 

as the equation of the real part, and 

2 2 

(2 u> g) X 2 + (- g - (o (g + (■£■) H g - (■£•) J) + A (^) + A (^-) g) X 

2 2 

+(£)(- A - A H {£) - J + E D + F B + F C (i) ) = 0 


(16) 


(17) 


( 18 ) 
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as the equation of the imaginary part. 

By solving these equations by the quadratic formula for X for assumed 
values of 1/K and plotting the square root of these roots against the 

parameter i = — , the frequency ratio (to /w) can be determined at 

NO) CX 

the intersection of the curves. 

Thus the condition for which flutter can occur is obtained from 
the relation 


a 

(19) 

^ X 


a) 


a 

(20) 


K~V X 


Solution for Flutter Mode Shape 

The mode shape for flutter is obtained by substituting uip and Vp 
into equation (2), normalizing on either q 1 or q 2 (say q 2 = 1.0) and 
solving for qj 



SOLUTION FOR DIVERGENCE SPEED 


The solution for the divergence speed is obtained by setting 
u = 0, g = 0 in equation (3) which gives 


V 11 


V 2 D 


12 


0 K 22 " V D 22 


= 0 


Expand (21) and solving for V gives 


K. 


V Div ~ 'j 


22 

} 22 


which in this case characterizes a torsional divergence problem. 

The divergence mode shape is obtained by substituting V div in 
equation (2) for w = 0, g = 0 and normalizing on say = 1*0 anc * 
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INPUT DESCRIPTION 


The "Namelist" method of input is used; see the Fortran Manual 
of Control Data System 6400/6600 Series for format instructions. 

The mode of all variables input is R-real. 

1. Namelist: $NAM1 including the following: 

Fortran name Description 

M Array defined as: 


DD 


R 




-ac 
a 2 + 


I 

m 


Aerodynamic matrix associated 
with displacements; defined as: 

r "i 


DD = | p c U L 

a 


0 

0 


c COS y 
b COS y 


Aerodynamic maxtrix associated 
with rates, defined 


R = p c Z C, 

a 




SEND NAM1 
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2. Namelist: $NAM2 
Fortran name 
G 

WALPH 


Wy 


Description 


Structural Damping, g 

a) = uncoupled natural frequency 

Oi 

defined as: 



to = uncoupled natural frequency 

V 

defined as: 


~\ f 


K (1, 1 

MUTT 


$ END NAM2 
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Program Flow Diagram 


C START ) 



• • SOLVE 
j FOR A, B, 

;c, D, E, F, H, J 


r X = 4.0 


. PRINT 

j HEADING 
• AND 
1 CASE NO. 


1 = 1 ! 


! XX (I) = X 
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SOLVE 

FOR 


lAl, Bl, Cl 



i 

l 

SOLVE 
FOR 
XI, X2 



SO 

F 

A2, B 

LVE 

OR 

2, C2 



S XI = XI j 

! 


S X2 = X2 

S X3 = X3 i 
1 


PRINT i 
U, S XI, S X2 j 
S X3 I 


R2 (I) = S X2 j 
R3 (I) = S X3 : 
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GP AND GN 

, > 


| XX (I) = X 

I 

; 1 — 


PRINT 

X, ZP, ZN, 
GP^ GN 

I 

G1 (I) = GP 

WW (I) = 
REAL (ZP) 
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( START ^ 
V FX J 


I SOLVE 
! FOR 
|A1, Bl, Cl j 


SOLVE 
FOR 
XI, X2 


SOLVE 

FOR 

A2, B2, C2 

, 



r 


. / 

\ 

G>0 

'nY 

N Y 


S 

II 

1 — 1 

X 

XI 

i 

S 

X2 = 

CM 

X 

S 

X3 = 

X3 



X3 = 

- B2 + (B2) 2 
2 A2 


4A2-C2 


! 
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( START ^ 
V FY J 

_ix: 

SOLVE FOR | 
AA , BB , CC | 


SOLVE 
FOR 
_ ZP_ 


i SOLVE 
i FOR 

! GP , 

I i 

i FY = GP - G 




SAMPLE PROBLEM 


The sample problem presented in this paper is that of a Pylon 
(strut) that was installed in the NASA Langley Research Center VSTOL 
tunnel for the purpose of housing wires running from a slip ring 
assembly through the tunnel wall to the recording instrumentation 
(see figure 1). 


Input Data 

The input data for this case is given by the listing on page 33. 

Output Data 

The output data consists of the tabulation of roots (both real 
and complex) for given values of 1/K (see pages 34 through 43). 
Following the roots tabulation the flutter frequency and flutter 
velocity are printed out. The program plot routine provides the 
parametric plots for determining the flutter condition from each of 
the two solution procedures. Plots are given in figures 2 through 
4 for a damping value of g = 0 and figures 5 through 7 for g = .01. 

The comparison of flutter speeds computed by the two methods 
and the root locus analysis for g = 0 is given as follows. 

Method 1 10872.9 

Method 2 10919.2 

Root Locus 10878.0 

As can be seen there is excellent agreement between the two methods as 

you would expect. These solutions are also compared to those obtained 

by performing a root locus analysis of the roots obtained by Program 

No. A4813 which also gave excellent agreement. 


“ F 

528.34 

528.34 

528.0 
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CONCLUDING REMARKS 

This paper has presented two methods for calculating flutter 
frequency and flutter speed for a system characterized by two degrees 
of freedom. These methods are simple to apply since the flutter speed 
and frequency can be solved for directly without resorting to a more 
complex eigenvalue approach such as a root locus analysis. Program 
A4868 is available through COSMIC. 
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APPENDIX A 

DEVELOPMENT OF GOVERNING EQUATIONS 
by 

William F. Hunter 

The equations of motion are derived for a simple two-degrees-of- 
freedom representation of the system. The bending and torsional 
deformations are expressed in terms of assumed displacement functions 
which satisfy the boundary conditions of the treated problem. The 
equations are developed for analyzing a strut in the VSTOL wind tunnel; 
but, with modification, they may be applied to an aircraft wing. 

The geometry of the swept strut, which has a uniform cross-section, 
is shown in the sketch below. 


z 



~ — center of pressure 
(a) Cross-sectional view (b) 



Sketch 1 
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The displacements z and a are assumed to vary with x according to 


| (x, t) = q 1 (t) sin 


a (x, t) = q 2 (t) sin 


( 1 ) 


such that q-| and q 2 become the generalized coordinates of the problem. 
The kinetic and strain energy are given by 


T = 




ca'Vfi + <H> 2 3 d * 

2 2 

[ El '(- — i— ) + G J (ff) ] d x 

3 XT 3 X 


( 2 ) 


where m is the distributed mass, I is the distributed mass moment 
of inertia about the center of mass axis, El is the bending stiffness, 
and GJ is the torsional stiffness. Noting that z = z - a a and 
letting m = m i and I = I cm i, the kinetic and strain energy may be 
expressed in quadratic form as 


T = \ (q} T [M] {q} 
U = \ (q) [K] {q} 


( 3 ) 
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(4) 


(5) 


( 6 ) 


Expressions for the generalized aerodynamic forces are obtained by 
considering a streamwise cross-section (see sketch 2) rather the 
previously depicted section which was taken normal to the x-axis. 



Sketch 2 
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Note that 

A- 

C = C/COS Y 

b = b/cos y 

(7) 

X = X COS Y 
a = a COS y 

Since the velocity of the center of pressure in the direction of 

z is 

A A 

Z„ = Z + b a 
Cp 

the total angle of attack is 

dj. = a - ^ C p/V - ol - z/V - b a /V (8) 

The aerodynamic loadings on the differential segment 





are approximated using simple strip theory by 



26 


L = 1 p V 2 S C L 
M = \ p V 2 S C M 


(9) 


where 

C L = C L “t 

a 

C M = C L b (10) 

S = c d x = cdx 


The virtual work of the aerodynamic forces acting on the differential 
segment is given by 

rw = L 6 + M <S a (11) 


substituting equations (1), (7), (8), (9), and (10) into equations 
(11) gives 


O' = f 


V 2 c C 


L cos Y - f q, - V q 2 > (c 6 q, + b S q 2 ) sin 2 (^) dx 


The virtual work for the entire strut is given by 


( 12 ) 


z 

6 VI = / 6 W 

0 

Integrating and comparing the result with 6 W = z Q. 6 defines the 

q i 

generalized forces Q-| and Q 2 - These forces may be expressed in a matrix 
equation as 



where the aerodynamics matrices [D] and [R] associated with the 
displacements and rates, respectively , are 


and 


[D] 


^ P C £ 


l 

V — 


c cos y! 

i 

b cos y! 

I 


[R] = ~ p c £ C L 

CL 



’i 


D C | 



_j 


{Ql 


j Q 1 


V 


( 14 ) 


(15) 


(16) 


Now that the kinetic energy, strain energy, and generalized forces 
have been expressed, Lagrange's equations yield 


[M] {q} + [K] {q} = {Q} (17) 

Introducing structural damping g and substituting for {Q } gives the 
governing differential equation as 
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PARTIAL PROGRAM LISTING 


PROGRAM FDETER ( INPUT , OUTPUT , TAPES , TAPE6 ) 
C PINNED-PINNED CASE 
REAL J,M(2.2>.KK<2,2) 

DIMENSION R<2,2), DD(2,2) 

DIMENSION R2<42>, R3<42). XX(42) 
DIMENSION EPS (3) 

REAL K 

COMPLEX ZP.ZN.BB.CC 
DIMENSION G1 ( 42 ) . UU(42) 

COMMON /BLK1/ A.B.C.D.E.F.H.J.U.G 
COMMON /BLK2/ RCROSS 
EXTERNAL FY 
EXTERNAL FX 

NAMELIST /NAM2/ G.UALPH.UY 
NAMELIST /NAM1/ M , DD , R 
DATA EPS/3* 001/ 

CALL PSEUDO 

ICASE-0 

READ (5/NAM1 ) 

IF ( EOF (5 ) ) 5.2 
READ (5.NAM2) 

ICASE- ICASE+1 
A-R( 1 , 1 )/M (1.1) 

B«M( 1 . 2 )/M( 1 . 1 ) 

C»DD( 1 . 2 )/M ( 1 . 1 ) 

D-RCl ,2)/M( 1 . 1 ) 

E-M<2. 1 )/M(2.2) 

F-R(2. 1 )/M( 2 . 2 ) 

H»DD(2.2)/MC2.2) 

J-R(2.2)/M(2.2) 

U-CUY/UALPH) 

X-4 0 

URITEC6. 11 )ICASE 
DO 3 1-1,40 
X-l/K 
AA«U**2 

S— (U**2*(l .-H*X**2)-1 . ) 
T--(A*X+U**2XJXX) 


A 1 
A 2 
A 3 
A 4 
A 5 
A 6 
A 7 
A S 
A 9 
A 10 
A 11 
A 12 
A 13 
A 14 
A 15 
A 16 
A 17 
A IS 
A 19 
A 20 
A 21 
A 22 
A 23 
A 24 
A 25 
A 26 
A 27 
A 28 
A 29 
A 30 
A 31 
A 32 
A 33 
A 34 
A 35 
A 36 
A 37 
A 38 



PARTIAL PROGRAM LISTING (CONT'D) 


BB»CNPLX(S.T) 

U- ( 1 . +HXX**2-A*JSXS*2-E*B-CS£*X**2+D*F*X**2 ) 
U»(X*(-A-ASH*X**2~J+E*D+B*F+C*FSX*S2) ) 

CC * CMP LX ( U . U ) 

ZP- (-BB+CSQRT(BB**2-4. *AA*CC ) )/<2 . *AA ) 

ZN- ( -BB-CSQRT ( BBS*2-4 . 2AA£CC ) )/ ( 2 . SAA ) 
GP-(AIMAG(ZP) )/(REAL(Z?>) 

Gff « (AIM AG ( ZN ) >✓ ( REAL ( 2,N ) ) 

XX(I )«X 

URITE ( 6 . 12 )X , ZP . ZN . GP. GN 

GKI )-GP 

UU( I )-REAL(ZP ) 

X-X+.5 

3 CONTINUE 

CALL RMULL (FY. 2 . , . 5. EPS, K . NUN, IERR ) 

X-K 

AA»U**2 

S«-(U**2*(1 ,-H*X**2)-l. ) 

T--(A*X+U**2*J*X) 

BB-CMPLX(S.T) 

U- ( 1 . +H*X**2-A*J*X**2-E*B-C*E2X2*2+D*F*X**2 ) 
U-(X*(-A-A*H*X**2-J+E*D+B*F+C*F*X**2) ) 
CC-CMPLXCU.U) 

ZP- <-BB+CSQRT(BB**2-4 *AA*CC ) )/(2 *AA ) 

GP- ( AIf1AG(ZP ) )/(REAL(ZP ) ) 
UU2-UALPH**2/(REAL(ZP)) 

UUU-SQRTC ABS(UU2 ) ) 

U«UUU*X 

URITE ( 6 , 13 )G . K 
URITE(6. 10)UUU.U 
CALL PLOT (XX. Gl, 1.2,1) 

CALL PLOT (UU.G1.3.2.2) 

X * 4 . 0 

URITE (6. 6 )ICASE 
DO 4 1-1.40 
XX(I )-X 
C X-l/K 

Al-U**2*(l.-G**2) 


A 39 
A 40 
A 41 
A 42 
A 43 
A 44 
A 45 
A 4S 
A 47 
A 43 
A 49 
A 50 
A 51 
A 52 
A 53 
A 54 
A 55 
A 56 
A 57 
A 58 
A 59 
A 60 
A 61 
A 62 
A 63 
A 64 
A 65 
A 66 
A 67 
A 68 
A 69 
A 70 
A 71 
A 72 
A 73 
A 74 
A 75 
A 76 


GO 

O 



o o o u) r- oo o> 


PARTIAL PROGRAM LISTING (CONT'D) 


4 


S' 


10 

11 


B1 ■ (-1 . -U**2*( 1 . +X**2*H+G*J*X ) ) 

Cl-1 . +Xt&2%H-Xtt3tAXJ-EtB-EtCtXtt2+FtDtXtX2 
Xl-<-Bl+SGRT(Bl**2-4 *A1*C1 ) )/(2. *A1 ) 
X2-(-Bl-SQRT(Bl**2-4.*Al*Cl ) )/C2.«Al ) 

A2-2. *U*G 

B2- ( -G-U**2* ( G+X**2*H*G-X*J )+A*X+A*X*G ) 
C2«X*(-A-A*H*X**2-J+E*D+F*B+F*C*X**2 ) 

X3--02/B2 

IF (G.GT.0.) X3- <-B2+SQRT(B2**2-4 . *A2*C2 ) )/(2 *A2 ) 
SX1 *SQRT(X1 ) 

SX2-SGRTCX2) 

SX3-SQRT(X3) 

URITE ( 6 . 7 )X , SXi , SX2, SX3 

R2( I )*SX2 

R3( I )*SX3 

X-X+.5 

CONTINUE 

CALL RMULL ( FX , 2 . , . 5 , EPS * < , NUM . IERR ) 

URITE(6,8)G,K 

CALL PL0T1 (R2.R3.XX) 

URITE ( 6 , 9 )RCROSS 
UUU-UALPH/RCROSS 
U-UUU*K 

UR I TEC 6, 10 )UUU,U 
GO TO 1 

CALL CALPLT (00,0.0,999.) 

STOP 


FORMAT ( 1H1 , 47H ROOTS OF REAL AND IMAGINARY EQUATIONS FOR CASE, 12) 
FORMAT ( / . 5H 1/K- , F6 . 3 , 3X , 7HSQRTX1 ■ , F15 . 6 , 3X , 7HSQRTX2- , FIS . 6 , 3X , 8H 
1 SQRTX3-.F1S.6) 

FORMAT (// , 3H G-,FS.2,9H FOR 1/K«,F8.4) 

FORMAT (/, 58H REAL AND IMAGINARY ROOTS INTERSECT UHEN SQUARE ROOT 
10F X-.F8.4) 

FORMAT (/, 3H U- , F12 . 4 , 2X . 2HU- , F12 . 4 ) 

FORMAT ( 1H1 , 4SH COMPLEX ROOTS OF FLUTTER DETERMINANT FOR CASE, 12) 


A 77 
A 78 
A 79 
A 80 
A 81 
A 82 
A 83 
A 84 
A 85 
A 86 
A 87 
A 88 
A 89 
A 90 
A 91 
A 92 
A 93 
A 94 
A 95 
A 96 
A 97 
A 98 
A 99 
A 100 
A 101 
A 102 
A 103 
A 104 
A 105 
A 106 
A 107 
A 108 
A 109 
A 110 


A 113 
A 114 



PARTIAL PROGRAM LISTING (CONCLUDED) 


12 

13 


a. 


FORMAT <//.SH 1/K" . F6 . 3.3X. 3HZP" . FJ5 . 6. IX. F15 . 6, 3X. 3HZN- . F15 . 6, IX. 
IF 15 . 6. 3X. 3HGP* . F15 . 6 . 3X. 3HGN* . F15 . 6 ) 

FORMAT <//.3H G-.F5.2.9H FOR 1/K-.F8.4) 

END 

FUNCTION FX (X) 

COMMON /BLK2/ RCROSS 
REAL J 

COMMON /BLK1/ A, B . C. D. E.F. H. J. U, G 
Al-U**2*(l .-G**2) 

B1 » < -1 . -U**2* ( 1 . +X**2*H+G*J*X ) ) 

Cl-1 . +X**2*H-X**2*A*J-E*B-E*C*X**2+F*D*X**2 
XI » ( -Bl+SQRT ( Bl**2-4 . *A1*C1 ) )/<2. *A1 ) 

X2» (-Bl-SGRT(Bl**2-4 . *A1*C1 ) )/C2 . *A1 ) 

A2-2.*U*G 

B2-(-G-U**2*(G+X**2*H*G-X*J)+A*X+A*X*G> 

C2«X*<-A-A*H*X**2-J+E*D+F*B+F*C*X**2) 

X3»-C2/B3 

IF (G.GT.0.) X3» < -B2+SQRT(B2**2-4 . *A2*C2 ) )/(2. *A2 ) 

SX1 “SORT ( XI ) 

SX2-SQRTCX2) 

RCR0SS-SX2 

SX3«SQRT(X3) 

FX-SX2-SX3 

RETURN 

END 

FUNCTION FV (X) 

COMMON /BLK1' A.B.C.D.E.F.H. J.U.G 
COMPLEX BB.CC.2P 
REAL J 
AA-U**2 

S--(U**2*<1 .-H*X**2>-1 . ) 

T— <A*X+U**2*J*X) 

BB-CMPLX(S.T) 

U- ( 1 . +H*X**2-A*J*X**2-E*B-C*E*X**2+D*F*X**2 ) 

U» (X* ( -A-A*H*X**2-J+E*D+B*F+C*F*X**2 ) ) 

CC-CMPLX(U.U) 

ZP- <-BB+CSQRT(BB**2-4 . *AA*CC ) )/(2 . *AA ) 

GP-(AIMAG(ZP))/(REAL(ZP)) 

FV-GP-G 

RETURN 


A 115 
A 116 
A 117 
A 118- 
B 1 
B 2 
B 3 
B 4 
B 5 
B 6 
B 7 
B 8 
B 9 
B 10 
B 11 
B 12 
B 13 
B 14 
B 15 
B 16 
B 17 
B 18 
B 19 
B 20 
B 21- 
E 1 
E 2 
E 3 
E 4 
E 5 
E 6 
E 7 
E 8 
E 9 
E 10 w 
E 11 N 
E 12 
E 13 
E 14 
E 15 



INPUT DATA FOR SAMPLE PROBLEM 
CASE 1 


.OOOsbi 

— 0 . * dKw' ) >1 s 3J 1 . * .i Y=£ i . J 7h 


INPUT DATA FOR SAMPLE PROBLEM 
CASE 2 

z,\^y : 1 >;sH.‘j02»-1.0i»-l.J'3«l«0 *0.* .00»l ♦0*00031 »^ a *ODft» .OOS6» .OliS6* 

• J 0 L* 7 o * 

i \ A M 2. i j = • 0 1 * * A l. <*i = v. u 1 • :> ^ 't ~ - 1 • f "• 



CO 

CO 



OUTPUT DATA FOR SAMPLE PROBLEM 


COMPLEX ROOTS OF FLUTTER DETERMINANT FOR CASE 1 


1/K- *.000 

ZP* 

-.927339 

.007338 

ZN* 

-790.382780 

1/K- 4.500 

ZP* 

-.931223 

.008168 

ZN * 

-790.380885 

1/K* 5.000 

ZP- 

-.935564 

.008968 

ZN* 

-790.378767 

1/K* 5.500 

ZP- 

-.940360 

.009733 

ZN* 

-790.376427 

1/K* 6.000 

ZP* 

-.945611 

.010461 

ZN* 

-790.373866 

1/K* 6.500 

ZP* 

-.951317 

.011149 

ZN* 

-790.371084 

i/K* 7.000 

ZP* 

-.957478 

.011792 

ZN* 

-790.368081 

1/K* 7.500 

ZP* 

-.964094 

.012388 

ZN* 

-790.364858 

1/K- 8.000 

ZP* 

-.971162 

.012933 

ZN* 

• 

-790.361416 

1/K* 8.500 

ZP* 

-.978684 

.013423 

ZN* 

-790.357754 

1/K- 9,000 

ZP* 

-.986659 

.013856 

ZN* 

-790.353874 

1/K* O’. 500 

ZP* 

-.9950755 • a 

.014228 

ZN- 

-790.349775 

l/K-10.000 

ZP* 

-1.003963 ~ 

' .014535 

ZN* 

-790.345459 

/ 

1/K* 10. 500 

ZP- 

. -1.013292 

.014775 

ZN* 

-790,340926 

U/K-11.000 

ZP* 

-1.023071 

.014944 

ZN= 

-790.336177 

1/K-.11.500 

ZP* 

-1.033299 

.015038 

ZN* 

-790.331213 


15.697760 

GP- 

-.007912 

GN* 

-.019861 

17.660067 

GP* 

-.008771 

GN- 

-.022344 ' 

19.622405 \ 

GP* 

-.009585 

GN- 

. -.024827 

21.584776 

GP* 

-.010351 

GN* 

-.027309 

23.547185 

GP* 

-.011063 

GN* 

-.029792 

25.509635 

GP* 

-.011719 

GN* 

-.032276 

27.472129 

GP* 

-.012316 

GN* 

-.034759 

29.434671 

GP* 

-.012849 

GN* 

-.037242 

31.397263 

GP* 

-.013317 

GN* 

-.039725 

33.359910 

GP* 

-.013715 

GN* 

-.042209 

35.322614 

GP* 

-.014043 

GN* 

-.044692 

37.285380 

GP- 

-.014298 

i 

GN* 

-.047176 

39.248209 

GP* 

-.014478 

GN* 

-.049660 

41.211107 

GP* 

-.014581 

GN* 

-.052143 

43.174075 

GP* 

-.014607 

GN* 

-.054627 

45.137118 

GP* 

-.014554 

GN* 

-.057112 

co 





-p» 


• 

1 

* 




OUTPUT DATA FOR SAMPLE PROBLEM (CONT'D) 


1/K*12 • OOO 

ZP* 

-1.043976 

.015055 

ZN = 

-790.326033 

47. 100239 

GP* 

-.014421 

GN* 

-.059596 

1/K* 12 • 500 

ZP* 

-1.055101 

. 014991 

ZN* 

-790.320640 

49.063440 

GP* 

-.014208 

GN* 

-.062080 

1/K* 13. 000 

ZP* 

-1.066673 

.014842 

ZN = 

-790.315034 ^ 

51.026726 

GP* 

-.013915 

'GN* 

-.064565 

1/K*13.500 

ZP* 

-1.078691 

.014606 

ZN * 

-790.309215 

52.990099 

GP* 

-.013541 

GN* 

-.067050 

1/K* 14. 000 

ZP* 

-1.091155 

.014280 

ZN* _ „ 

'-790. 303185 

54.953563 •. 

GP* 

-.013087 

GN* 

* -.069535 

i/K*14.500 

ZP* 

-1.104063 

.013859 

ZN* 

-790.296944 

56.917121 

GP* 

-.012553 

GN* 

-.072020 

1/K*15.000 

ZP* 

-1.117415 

.013341 

ZN* 

-790.290494 

58.880776 . 

GP* 

-.011939 

GN* 

-.074505 

1/K*15. 500 

ZP* 

-1.131209 

.012722 

ZN* . 

-790.283835 

60.844532 

GP* 

-.011246 

* GN* 

-.076991 

1/K*I 6 .000 

ZP* 

-1.145445 

.011999 

ZN* 

-790.276968 

62.808392 

GP* 

-.010476 

GN* 

-.079476 

1/K*16.500 

ZP* 

-1.160122 

.011170 

ZN* 

-790.269894 

64.772359 

GP* 

-.009628 

GN* 

-.081962 

1/K=17.000 

ZP* 

-1. 175233 

.010230 

ZN* * 

-790.262615 

66.736436 

GP* 

-.008705 

GN* 

-.084448 

1/K* 17. 500 

ZP* 

-1.190793 

.009177 

ZN* 

-790.255132 

68.700626 

GP* 

-.007706 

GN* 

-.086935 

1/K*1 8. 000 

ZP* 

-1.206785 

.008007 

ZN* 

-790.247444 

70.664934 

GP* 

-.006635 

GN* 

-.089421 

i/K*18.500 

ZP* 

• 

-1.223214 

.006717 

ZN* 

-790.239555 

72.629360 

GP* 

-.005491 

i 

GN* 

-.091908 

1/K*1 9. 000 

ZP* 

-1.240077 \ 

.005304 

ZN* 

-:790. 231464 

74.593910 

GP* 

-.004278 

GN* 

-.094395 

1/K*19. 500 

ZP* 

-1.257375 

.003766 

ZN* 

i 

-790.223174 

76.558586 

GP* 

. -.002995 

GN* 

-.096882 

1/K*20. 000 

ZP* 

. -1.275105 

.002098 

ZN* 

-790.214684 

78.523392 

GP* 

-.001645 

GN* 

-.0993 70 


A 


CO 

CJT 



OUTPUT DATA FOR SAMPLE PROBLEM (CONT'D) 


1/K = 20.500 

ZP- 

-1# 293266 

.000297 

ZN = 

-790.205998 

00.488329 

GP-- 

* -.000230 

GN= 

-.101857 

1/K- 21 • 000 

ZP- * 

. -1.311858 

-.001639 

ZN = 

-790.197115 

82.453402 

GP- 

.001249 

’ GN- 

-.104345 

I /K*21 *500 

ZP- 

-1. 330878 

-.003713 

ZN = 

-790.188037 

84.418614 

GP- 

.002790 

GN- 

-.106834 

1/K*22.000‘ 

ZP- 

-1.350326 

-.005929 

ZN= 

-790.178766 

86.383967 

GP* 

.004391 

GN- 

-.109322 

l/K-22 • 500 

ZP- 

-1.370200 

-.008290 

ZN = 

-790. 169302 

88.349465 

GP- 

.006050 

GN- 

-.111811 

1/K*23. 000 

ZP* 

• 

-1.390499 

-.010798 

ZN a 

-790.159648 

90.315111 

GP- 

.007766 

GN- 

-.114300 

1/K=23. 5 00 

ZP- 

-1.411221 

-.013458 

ZN = 

-790.149804 

92.280908 

GP- 

.009536 

GN- 

-.116789 


G= 0*03 FOR 1/K- 20.5791 
W* 528.3450 V- 10872.8658 

I 




OUTPUT SAMPLE PROBLEM 

ROOTS OF REAL AND IMAGINARY EQUATIONS FOR CASE 1 


1/K* 4.000 

SQRTX 1= 

23.149339 

SQRTX2 » 

.961691 

SQRTX3* 

1.138499 

1/K* 4.500 

SQRTX 1 = 

28.149306 

SQRTX2 = 

.963635 

SQRTX3* 

1.138499 

' 1/K* 5.000 

SQRTX 1 3 

28.149269 

SQRTX2= 

.965907 

SQRTX3* 

1. 138499 

1/K* 5.500 

SQRTX 1* 

28.149229 

SQRTX2 = 

.968358 

SQRTX3* 

1.138499 

1/K* 6.003 

SQRTX 1* 

28.149184 

SQRTX2= 

.971036 

SQRTX3* 

1.138499 

l/K* 6.500 

SQPTX1- 

28.149136 

SQRTX2= 

.973937 

SQRTX3* 

1.138499 

l/K« 7.000 

SQRTX1* 

23.149034 

SQRTX2 3 

.977062 

SQRTX3* 

1.138499 

1/K* 7.500 

SQRTX 1 * 

28.149028 

SQRTX2* 

.980406 

SQRTX3* 

1.133499 

1/K* 8.000 

SQRTX 1 * 

28.148968 

SQRTX2” 

.983969 

SQRTX3* 

1.138499 

1/K* S'. 503 

SQRTX 1* 

' 28.148904 

SQRTX2 = 

.987747 

SQRTX3* 

1.138499 

1/K= 9.003 

SQRTX la 

28.148836 

SQRTX2 3 

.991739 

SQRTX3* 

1.138499 

l/K= 9.500 

SQRTX 1 * 

28.143765 

SQRTX2= 

.995942 

SQRTX3* 

1.138499 

1/K= 10.000 

. SQRTX 1» 

23.143689 

SQRTX2* 

1.000352 

SQRTX3* 

1.138499 

1/K=1 0. 500 

SQRTX 1* 

28.143610 

SQR TX2- 

1.004968 

SQRTX3* 

1.138499 

1/K=1 1 .000 

SQRTX 1* 

28.148527 

SQRTX2 3 

1.009787 

SQRTX3* 

1.138499 

1/K* 11 . 503 

SQRTX 1* 

28.148440 

SQRTX2= 

1.014805 

SQRTX3* 

1.138499 

1/K=12. 000 

SQRTX 1 * 

28.148349 

SQRTX2- 

1.020020 

SQRTX3* 

1.138499 

l/K-12.500 

SQRTX 1* 

23.143255 

SQRTX2® 

1.025429 

SQRTX3* 

1.138499 

1/K*13. 000 

SQRTX 1= 

28.148156 

S0RTX2= 

1.031028 

SQRTX3* 

1.138499 

1/K*13. 500 

SQRTX 1* 

28.148054 

S0RTX2* 

1.036615 

SQRTX3* 

1.138499 

1/K*! 4.000 

SQRTX 1» 

28.147947 ' 

SQRTX2* 

1.042787 

SQRTX3* 

; 

1.138499 

1/K* 14. 500 

SQRTX 1- 

28.147837 

$QRTX2= 

1.048940 

SQRTX3® 

1.138499 

i/K*15 .000 

SQRTX 1* 

23.147723 

SQRTX2=* 

1.055271 

SQRTX3* 

• 1.138499 

1/K*15. 500 

SQRTX 1 * 

28.147605 

SQRTX2* 

1.061777 

SQRTX3- 

1.138499 



1/K= 16.000 

SQRTX 1 = 

28.147483 

SQRTX2= 

1/K=: 6.500 

SOR T X 1 = 

28.147358 

SQRTX2= 

1/K= 17. 000 

S0R T X 1= 

2 8.147228 

SQRTX2= 

1/K=) 7. 500 

SpRi-yi* 

28.147095 

SQRTX2= 

1/K= 13.000 

SQRTX1* 

23.146957 

SQRTX2= 

1/K= 18.500 

SQRTX1= 

28.146316 

SQRTX2= 

1/K* 19. 000 

SQRTX is 

28.146671 

SQRTX2 3 

1/K*19. 500 

SQRTX 1 * 

28.146522 

SQR TX2= 

1/K=20. 000 

SQRTX 1= 

28.146370 

SQRTX2= 

1/K=20. 500 

SQRTX 1= 

23.146213 

SQRTX2= 

1/K=21 .000 

SQRTX 1= 

29.146052 

SQRTX 2 = 

1/K = 21. 500 

SQRTX 1 = 

20.145888 

SQRTX 2= 

1/K=22. 000 

SQRTX 1 = 

28.145720 

SQRTX2= 

1/K=22.500 

SQRTX 1= 

28.145547 

SQRTX2* 

l/K=23.900 

SQRTX 1* 

28.145371 

SQRTX2= 

1/K=23. 500 

SQRTX 1= 

23.145191 

SQRTX2= 

G= 0.00 FOR 

1 At- 20.6671 




RE4L AND IMAGINARY ROOTS INTERSECT WHEN SQUARE 
W* 528.3405 V« 10919.2489 


OUTPUT SAMPLE PROBLEM (CONT 

•D) 

1.069455 

SQRTX3 = 

1.138499 

.1.075302 

SQRTX3= 

1.138499 

1.082314 

SQRTX3= 

1. 138499 

1.089489 

SQRTX3= 

1.138499 

1.096822 

SQRTX3= 

1.138499 

1. 104312 

SQRTX3= 

1.138499 

1. 111954 

SQRTX3* 

1. 138499 

1.119746 

SQRTX3 = 

1. 138499 

1.127684 

SQRTX3 * 

1.138499 

1.135767 

SQRTX3= 

1.138499 

1.143999 

SQRTX3 = 

1.138499 

1.152349 

SQRTX3= 

1.138499 

1. 160844 

SQRTX3= 

1.138499 

1.169471 

SQRTX3 s 

1.138499 

1.178226 

SQRTX3* 

1.138499 

1. 187107 

SQRTX3= 

1.138499 


ROOT OF X = 1.1385 


OUTPUT DATA FOR SAMPLE PROBLEM 


COMPLEX ROOTS OF FLUTTER DETERMINANT FOR CASE 2 


1/K* 4.000 

ZP* 

-• 927339 

. .007338 

ZN* 

-790.382780 

15.697760 

GP- 

-.007912 

GN* 

-.019861 

l/K* 4.500 

ZP* 

-•931223 

•008168 

ZN* 

-790.380885 

17.660067 

GP= 

-.008771 

GN* 

-.022344 

l/K* 5.000 

ZP* 

-.935564 

» 

.008968 

ZN* 

-790.378767 

19.622405 

GP* 

-.009585 

GN* 

-.024827 

• 

l/K- 5.500 

ZP* 

-•940360 

.009733 

ZN* 

-790.376427 

21.584776 

GP* 

-.010351 

GN* 

-.0273 09 

l/K* 6.000 

ZP- 

-•945611 

•010461 

ZN* 

-790.373866 

23.547185 

GP* 

-.011063 

GN* 

-.029792 

l/K* 6.500 

ZP* 

-•951317 

.011149 

ZN* 

-790.371084 

25.509635 

GP* 

-.011719 

GN* 

-.032276 

l/K* 7.000 

ZP* 

‘ -.957478 

.011792 

ZN* 

-790.368081 

27.472129 

GP* 

-.012316 

GN* 

-.034759 

l/K* 7.500 

ZP* 

-•964094 

.012388 

ZN* 

-790.364858 

29.434671 

GP* 

-.012849 

GN* 

-.037242 

l/K* 8.000 

ZP* 

-.971162 

.012933 

ZN* 

-790.361416 

31.397263 

GP* 

-.013317 

GN* 

-.039725 

l/K* 8.500 

ZP* 

-•978684 

.013423 

ZN* 

-790.357754 

*33.359910 

GP* 

-.013715 

GN* 

-.042209 

l/K* 9.000 

ZP- 

-.986659 

.013856 

ZN* 

-790.353874 

35.322614 

GP* 

• 

-.014043 

GN* 

-.044692 

4 

i/K* 9.500 
• 

ZP* 

-.995085 

*7 . . 

.014228 

ZN* 

-790.349775 

37.285380 

GP* 

-.014298 

GN* 

-.047176 

i/K*to;ooo 

'ZP- 

.1 * 

-1.0O3$63 

«. ■ 

.014535 

ZN* 

-790.345459 

39.248209 

GP* 

-.014478 

GN* 

-.049660 

1/K*1 0. 500 

ZP* 

-1.013292 

.014775 

ZN* 

-790.340926 

41.211107 

GP* 

-.014581 

GN* 

-.052143 

l/K* 11. 000 

ZP* 

-1.023071 

.014944 

ZN* 

-790.336177 

43.174075 

GP* 

-.014607 

GN* 

-.054627 

1/K*U,500 

ZP* 

-1,033299 

.015038 

ZN* 

-790.331213 

45.137118 

GP- 

-.014554 

GN- 

-.057112 



OUTPUT DATA FOR SAMPLE PROBLEM (CONT'D) 


1/K* 1 2. 000 

ZP* 

-1.043976 

. 015055 

ZN = 

-790.326033 

47. 100239 

GP* 

-.014421 

GN* 

-.059596 

l/<= 1 2. 5 00 


-1.055101 

. 014991 

ZN* 

— 790 . 32 0640 

49.063440 

. GP* 

-.014208 

GN* 

-.062080 

1/K= 13.000 

ZP* 

-1.066673 

. 014842 

ZN = 

-790 .31 5034 

51.026726 

GP* 

-.013915 

GN* 

-.064565 

1/K*13.500 

ZP* 

-1.078691 

. 014606 

ZN = 

-790.309215 

5*2. 990099 

GP* 

-.013541 

GN* 

-.067050 

1/K=14.000 

ZP* 

-1.091155 

.014280 

ZN* * 

-790. 303185 

54.953563 

GP* 

-.013087 

GN* 

-.069535 

1/K* 14*500 

ZP* 

-1. 104063 

.013359 

ZN = 

-790.2969*4 

56.917121 

GP* 

-.012553 

GN* 

-.072020 

1/K*15.000 

ZP* 

-1.117415 

.013341 

ZN* 

-790.290494 

58.880776 

GP* 

-.011939 

GN* 

-.074505 

1/K* 15.500 

ZP* 

-1.131209 

.012722 

ZN* 

-790.283835 

60.844532 

GP* 

-.011246 

GN* 

-.076991 

1/K*1 6. 000 

ZP* 

-1.145445 

.011999 

ZN* 

-790.276968 

62.808392 * 

GP* 

-.010476 

GN* 

-.079476 

1/K* 1 6* 500 

ZP* 

-1. 160122 

.011170 

ZN* 

-790.269894 

64.772359 

GP* 

-.009628 

GN* 

-.081962 

1/K*1 7. 000 

ZP* 

-1. 175239 

.010230 

ZN* 

-790.262615 

66.736436 

GP* 

-.008705 

GN* 

-.084448 

1/K=17.500 

ZP* 

-1.190793 

.009177 

ZN* 

-790.255132 

68.700626 

GP* 

-.007706 

GN- 

-.086935 

1/K=18.000 

ZP* 

-1.206785 

.008007 

ZN* 

-790.247444 

70.664934 

GP* 

„ -.006635 

GN= 

-.089421 

1/K*1 8* 500 

ZP* 

-1.223214 

.006717 

ZN* 

-790.239555 

72.629360 

GP* 

-.005491 

GN* 

-.091908 

1/K*1 9. 000 

ZP* 

-1.240077 

.005304 

ZN* 

-790.23^464 

74.593910 

GP* 

-.004278 

GN* 

-.094395 

i/K*19. 5 00 

ZP* 

-1.257375 

.003766 

ZN* 

-790.223174 

76.558586 

GP* 

-.002995 

GN* 

-.096882 

1/K*20.000 

ZP* _ 

-1.275105 

.002098 

ZN* 

-790.214684 

78.523392 

GP* 

-.001645 

GN* 

-.099370 


4 ^ 

O 



1 /K*20. 500 

ZP* 

-1.293266 

OUTPUT 

.000297 

DATA 

ZN = 

FOR SAMPLE PROBLEM 

-790.205998 

(CONT'D) 

80.488329 

GP* 

-.000230 

GN* 

1/K=21 .000 

ZP* 

-1.311958 

-.001639 

ZN 3 

-790.197115 

82.453402 

GP = 

.001249 

GN* 

1/K=21. 500 

ZP* 

-1.330878 

-.003713 

ZN= 

-790. 188037 

84.418614 

GP* 

.002790 

GN* 

l/K=22.000 

ZP* 

-1.350 326 

-.005929 

ZN* 

-790.178766 

86.383967 

GP* 

.004391 

GN* 

l/K-22.500 

ZP* 

-1.370200 

-.008290 

ZN = 

-790. 169302 

88.349465 

GP* 

.006050 

GN* 

1/K* 23.000 

ZP* 

-1.390499 

-.010798 

ZN* 

-790.159648 

90.315111 

GP* 

.007766 

GN* 

1/K*23.500 

ZP* 

-1.411221 

-.013458 

ZN* 

-790.149804 

92.280908 

GP* 

.009536 

GN* 


6* .31 FOR 1/K* 23.6285 

505.3868 V* 11941.5551 


-.101857 
-.1043 45 
-.106834 
-.109322 
-.111811 
-.114300 
-.116789 



OUTPUT DATA FOR SAMPLE PROBLEM (CONT'D) 

ROOTS OF REAL AMO IMAGINARY EQUATIONS FOR CASE 2 


1/K* A. OOO 

$0R T X1= 

23. 1 50/49 

SQR TX2 = 

.961691 

SQRTX3 = 

1.489313 

1/K* 4.500 • 

SQR TX 1 3 

28.150715 

SQR TX2 = 

.963684 

SQRTX3= 

1.437943 

1/K* 5.000 

SORTXl* 

28.150679 

S0RTX2* 

.965907 

$0RTX3= 

1.399236 

1/K* 5.500 

SQRTXl 3 

28.150638 

SQRTX2 s 

.968358 

SQRTX3= 

1.369109 

1/K= 6.000 

SQRTX1* 

28.150594 

SQRTX2= 

.971036 

SQRTX3= 

1.345031 

I/K= 6.500 

SQRTX l 3 

28.150546 

SQRTX2 3 

. 973937 

SQRTX3 s 

1.325365 

1/K= 7.000 

SQRTX1 3 

23.150494 

SQRTX2 = 

.977062 

SQRTX3 = 

1.309012 

1/K* 7.500 

SQP T X 1* 

28.15043 8 

SQRTX2= 

.900406 

SQRTX3= 

1.295204 

!/K= a. 000 

SQRTX1= 

28.150378 

SQRTX2= 

.903969 

SQRTX3- 

1.283393 

1/K= 8.500 

SQRTX1 3 

28.150314 

SQRTX2- 

.987747 

SQRTX3= 

1.273177 

1/K= 9.000 

SQRTX 1* 

2 8.150246 

SQRTX2= 

.991739 

SQRTX3 35 

1.264256 

1/K= 9.500 

SQtvTvC^ 

28.150175 

SQRTX2 = 

.995942 

SQRTX3= 

1.256397 

1/K = 10. 000 ' 

SQRTX1* . 

28.150100 

SQRTX2- 

1.000352 

SQRTX3= 

1.249424 

1/K=10. 500 

SQRTX 1 = 

28.150020 

SQRTX2* 

1.004968 

SQRTX3 3 

1.243194 

1/K=ll; 000 

SQRTX1 3 

28.149937 

SQRTX2= 

1.009787 

SQRTX3 s 

1.237595 

1/K=I1.500 

SQRTX 1= 

28.149851 

SQRTX2= 

1.014005 

SQRTX3= 

1.232535 

1/K= 12. 000 

SQRTX l 3 

28.149760 

. SQR TX2 = 

1.020020 

SQRTX3= 

1.227942 

I/K*12. 500 

SQRTX 1 “ 

28. 149665 

SQRTX2= 

1.025429 

SQRTX3 * 

1.223752 

1/K*13.000 . 

SQRTX 1 3 

28.149567 

SQRTX2- 

1.031028 

SQRTX3= 

1.219916 

I /K= 13. 5 00 

SQRTX 1= 

28.149464 

SQRTX2* 

1.036815 

SQRTX3 3 

1.216390 

1/K=14.000 

SQR T X I 3 

28.149358 

SQRTX2* 

1.042787 

SQRTX3- 

1.213138 

1/K*14.500 

SQRTX 1« 

28.149248 

SQRTX2 3 

1.048940 

SQRTX3® 

1.210130 

1/K=15. 000 

SQRTX1* 

28.149134 

SQRTX2= 

1.055271 

SQRTX3 3 

1.207339 

1/^*15.500 

SQRTX1* 

28.149016 

SQRTX2* 

1.061777 

SQRTX3* 

1.204742 


1 /K= 16. 000 

SQP T Xi= 

28. 148894 

OUTPUT 

SQRTX2 S 

DAT fto F ^ 5 i AMPl l ) «W. LEH 

1/K* 1 6. 500 

SQR T X 1= 

28.148769 

SQRTX2 S 

1.075302 

S QRTX3 = 

1.200057 

1/K=1 7 .000 

SQRTX 1* 

28.148639 

SQRTX2= 

1. 00231 A 

SQRTX3 = 

1.197936 

1/K=17.500 

SQRTX 1* 

28. 148506 

SQRTX2= 

1.089488 

SQRTX3 = 

1.195944 

l/K= 19.000 

SQRTX 1= 

28.149369 

S0RTX2= 

1.096822 

5QRTX3= 

1.194071 

1/K= 18. 500 

SQRTX1= 

28.143223 

SQRTX2= 

1.104311 ‘ 

SQRTX3= 

1.192306 

l/K= 19.000 

5QRTX1= 

28.148083 

SQPTX2= 

1. 111954 

SQRTX3 = 

1.190640 

1/K=1 9. 500 

SQRTX 1= 

28.147934 

SQRTX2= 

1. 119746 

SQRTX3= 

1.189065 

1/K=20. 000 

SQRTX 1= 

28.147781 

SQRTX 2= 

1.127684 

S0RTX3 = 

1.187573 

l/K=20. 500 

SQRTXl* 

28.147624 

SQRTX2= 

1.135766 

SQRTX3- 

1.186159 

1/K=2t. 000 

SQRTX1* 

• 28.14T464 

SQRTX2= 

1.143989 

SQRTX3 = 

1. 184815 

1/K*21.500 

SQPTX1* 

28.147300 

SQRTX2 3 

1. 152349 

SQRTX3= 

1.183538 

1/K= 22 • 000 

SQRTX 1= 

28.147131 

SQRTX2= 

1. 160844 

S0RTX3= 

1. 182322 

1/K=22. 500 

SQRTX 1= 

28.146959 

SQRTX2= 

1.169470 

SQRTX3* 

1.181163 

1/K=23. 000 

SQRTX 1* 

28.146783 

SQRTX2 3 

1.178225 

SQRTX3* 

1.180058 

1/K=23. 500 

SQRTX1= 

28.146603 

SQRTX2= 

1.187107 

SQRTX3® 

1.179001 


G* .01 FOR 1/K* 23.0925 

REAL AND IMAGINARY POOTS INTERSECT WHEN SQUARE ROOT OF X* 1.1799 
W* 509.9030 V* 11772. 6144 
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FIGURE I. Pylon Installation 
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FIGURE 5. 


Variation of 1/K with g 


- Solution Method (1) for Actual 
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FIGURE 6. Variation of Frequency Ratio with g - Solution Method (1) for Actual g.p .01. 
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